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The effectiveness of unsteady zero-net-mass-flux jets for fuselage drag reduction was evaluated nu- 
merically on a generic rotorcraft fuselage in forward flight with a rotor. Previous efforts have shown 
significant fuselage drag reduction using flow control for an isolated fuselage by experiment and nu- 
merical simulation. This work will evaluate a flow control strategy, that was originally developed on 
an isolated fuselage, in a more relevant environment that includes the effects of a rotor. Evaluation of 
different slot heights and jet velocity ratios were performed. Direct comparisons between an isolated 
fuselage and rotor/fuselage simulations were made showing similar flow control performance at a -3° 
fuselage angle-of-attack condition. However, this was not the case for a -5° angle-of-attack condition 
where the performance between the isolated fuselage and rotor/fuselage were different. The fuselage 
flow control resulted in a 17 % drag reduction for a peak Cp of 0.0069 in a forward flight simulation 
where p = 0.35 and CT /a = 0.08. The CFD flow control results also predicted a favorable 22% reduc- 
tion of the fuselage download at this same condition, which can have beneficial compounding effects on 
the overall performance of the vehicle. This numerical investigation was performed in order to provide 
guidance for a future 1/3 scale wind tunnel experiment to be performed at the NASA 14-by 22-Foot 
Subsonic Tunnel. 


Nomenclature 


Acs 

— 

fuselage body cross-sectional area 
(maximum), inches 

Aj 

= 

jet slot area, inches 

a<x> 

= 

freestream speed of sound, ft/s 

c D 

= 

fuselage drag coefficient, D / ( qooAcs ) 

C L 

= 

fuselage lift coefficient, L ( qooAcs ) 

C p 

= 

pressure coefficient, (p — p s ) / qoo 

Cp 

= 

rotor thrust coefficient, T / poo7tR 2 (£lR) 2 

Ct /g 

= 

blade loading coefficient 

Cp 

= 

jet momentum coefficient, YSjhjUj) / {q^Acs) 

D 

= 

drag force, lbs 

f 

= 

excitation frequency, Hz 

F + 

= 

reduced frequency, fW /Uoo 

h 

= 

jet slot height, inches 

H 

= 

fuselage main body height (maximum) 

L 

= 

lift force, lbs 

M 

= 

Mach number, Uoo/aoo 

P 

= 

pressure, psi 

Ps 

= 

static pressure, psi 

qoo 

= 

freestream dynamic pressure, 1 /2poo U 2 , psi 

R 

= 

rotor radius, inches 

R f 

= 

reference rotor radius, inches 

Re 

= 

Reynolds number, Uoo(2Rf)/v 

Uoo 

= 

freestream velocity, ft/s 
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U tip 

Uj 

VR 

w 

W 

x/Rp 

y/ R F 

z/R f 

ZNMF 

aj 

V 
a 

a j 

a s 

P 

V 

P 

Pi 

a 


rotor tip hover velocity, ft/s 
jet exit velocity, ft/s 
jet velocity ratio, Uj/Uoo 

width of surface blowing boundary condition, inches 

fuselage body width (maximum), inches 

normalized streamwise coordinate 

normalized cross-stream coordinate 

normalized vertical coordinate 

Zero-Net-Mass-Flux 

jet inclination angle, degrees 

jet phase angle, degrees 

rotor azimuth angle, degrees 

angle of attack, degrees 

jet inclination angle, degrees 

rotor shaft angle, degrees 

rotor advance ratio, Uoo cos a/ ( £IR ) 

kinematic viscosity, ft 2 /s 

density, slugs/ ft 3 

jet density, slugs/ ft 3 

thrust weighted rotor solidity 


Subscripts 

j = jet 

°o = freestream 


T he design of a rotorcraft fuselage shape places greater 
emphasis on function than on aerodynamic efficiency, 
typically resulting in designs that experience significant 
drag at cruise. The cruise drag of a rotary-wing aircraft is 



an order of magnitude higher than the cruise drag of a fixed- 
wing aircraft of the same gross weight (Ref. 1). At high ad- 
vance ratios, half the power to the main rotor is used to over- 
come the fuselage drag (Ref. 2). The drag force acting on 
a rotorcraft fuselage that features an aft- facing ramp is typ- 
ically dominated by massive flow separation on the rear of 
the fuselage. This results in a large pressure drag contribu- 
tion to the total drag force. The shear layers emanating from 
the separated region tend to roll-up in an unsteady fashion. 
This, in turn, results in unsteady aerodynamic loadings on 
the tail boom, stabilizers and tail rotor. By applying ac- 
tive flow control (AFC) to the ramp region of the fuselage, 
the size of the fuselage wake can be reduced, lowering the 
drag and reducing the undesirable unsteady loading on the 
downstream structures of the rotorcraft. Previous research 
in this area has demonstrated the importance of considering 
the impact of any drag reduction on the download that the 
fuselage develops. The works Martin, et al. (Ref. 3) and 
Ben-Hamou et al. (Ref. 4) demonstrated a fuselage drag re- 
duction through the application of active flow control, by 
20% and 11% respectively. Martin, et al. went on to show 
the profound effect that AFC can have on download reduc- 
tion. Since in a rotorcraft system the rotor is responsible for 
both thrust and lift, a benefit from a reduction in drag may 
be offset by an increase in download, effectually limiting 
any benefit from applying the flow control. 

In order to investigate different fuselage drag reduction 
approaches, an integrated experimental and computational 
research effort was undertaken under the NASA Subsonic 
Rotary Wing (SRW) Project utilizing a generic fuselage as 
its basis. A new generic fuselage was developed for this 
effort called the ROBIN-mod7. It is a modified version of 
the ROBIN (ROtor Body INteraction) fuselage that was de- 
veloped at NASA Langley in the 1970’s to be representa- 
tive of a generic helicopter and also to be easily reproduced 
for calculations (Ref. 5). The original ROBIN fuselage is 
well-documented in wind tunnel investigations (Refs. 6, 7) 
and is a popular generic fuselage within the rotorcraft CFD 
community (Refs. 8,9). The fuselage calculation procedure 
and the modified coefficients for the ROBIN-mod7 are dis- 
cussed fully by Schaeffler et al. (Ref. 10). It is intended that 
the ROBIN-mod7 fuselage can be established as a standard 
test fuselage for evaluating different AFC strategies both 
experimentally and computationally. 

The integrated experimental and computational ap- 
proach of this work allowed the baseline aerodynamic 
characteristics to be documented and validated simultane- 
ously (Ref. 10). Out of this effort and a related collaboration 
between NASA and ONERA, a flow control strategy was 
identified for the isolated fuselage which produced a 26% 
drag reduction using unsteady zero-net-mass-flux (ZNMF) 
jets (Ref. 10). The next step in this work is to see how 
the developed strategy is affected by increases in Reynolds 
number and by the effect of the rotor down wash. To ac- 
complish this, a wind tunnel test is planned in the NASA 


14-by 22-Foot Subsonic Tunnel. The research reported in 
this paper involves numerical simulations using a Reynolds- 
averaged Navier- Stokes (RANS) flow solver in support of 
this future wind tunnel entry. The investigation takes the 
lessons learned from the flow control strategies developed 
on the isolated fuselage and applies them to the more rele- 
vant environment of a fuselage and rotor operating at a high 
advance ratio. This allows for the study of the affect of 
the rotor and its downwash on the fuselage drag and down- 
load and on the effectiveness of the active flow control for 
elevating some of the drag and download forces. Addition- 
ally, results from this investigation will provide guidance on 
the flow control approach and identify slot configurations 
and flow control actuator parameters such as slot size, fre- 
quency, and peak jet velocity as well as to identify critical 
measurements needed for CFD validation. 

Numerical Modeling Approach 

Flow Solver 

The flow field for the fuselage was computed using 
the flow solver code, OVERFLOW (Refs. 11, 12) de- 
veloped at NASA. This code solves the compressible 
RANS equations using the diagonal scheme of Pulliam and 
Chaussee (Ref. 13). The RANS equations are solved on 
structured grids using the overset grid framework of Ste- 
ger, Dougherty, and Benek (Ref. 14). The flow is computed 
on node-centered grids composed of curvilinear body-fitted 
grids overset onto automatically generated Cartesian block 
background grids. Convergence and accuracy of the flow 
solution was improved using a Low-Mach number Precon- 
ditioning (LMP) option for steady simulations. The numer- 
ical simulations were performed using the parallel version 
of the OVERFLOW code developed by Buning (Ref. 15). 
This code uses the Message-Passing Interface (MPI) and 
can run on a tightly-coupled parallel machine or a network 
of workstations. The code distributes zones to individual 
processors and can split larger individual zones across mul- 
tiple processors using a domain decomposition approach. 

The RANS equations are solved implicitly using the 
Beam- Warming block tridiagonal scheme with a 3 nd order 
Roe upwind scheme for the inviscid flux terms. The Spalart 
Allmaras (SA) turbulence model was used for the numer- 
ical simulations (Refs. 16, 17). A Rotational/Curvature 
Correction (RCC) model, as implemented in OVERFLOW, 
was used to correct the SA model for surface curvature ef- 
fects (Ref. 18). The RCC model was shown in the previous 
effort to improve the symmetry of the baseline flow simula- 
tions and had a significant performance impact on simula- 
tions with small AFC slots. 

The blowing jets in OVERFLOW were simulated using 
a surface boundary condition, simplifying placement of the 
jets on the fuselage. Since the blowing slots intersect the 
fuselage surface at an angle, a 7 , the width of the jet surface 



boundary condition is: 

w = h/sin(aj) (1) 

where w is the width of the blowing boundary condition 
on the fuselage surface and h is the idealized flow control 
slot height. The jet blowing angle is based on the local sur- 
face tangent, normal to the slot span. An a j = 90° is there- 
fore normal to the local surface and a j = 0° is tangent to 
the local surface and normal to the slot local spanwise tan- 
gent. The jet properties were defined by setting a mass flow 
ratio per unit area of, PjUj/pooUoo, and the total tempera- 
ture ratio, T t j/T t oo. Unsteady ZNMF jets were simulated 
by modifying the surface boundary condition using a sinu- 
soidal fluctuation term defined as: 

Uj (t)=VR Uoo sin(27T ft + 4>) (2) 

where / is the frequency, (]) the phase, and VR is the peak 
jet velocity ratio. 

Grids 

The grids, shown in Figs. 1 and 2, included five fuselage 
grids with eight refined overset grids to resolve the flow 
control blowing on the surface and twelve grids for the four 
rotor blades. The rotor and body wakes were captured us- 
ing grids that have a uniform grid spacing equal to 10% 
of the rotor tip chord in all three directions. These wake 
grids extended 0.87? downstream of the rotor tip path plane 
and 0.87? below the rotor hub center, capturing the rotor and 
fuselage wake. Automatic background grids were gener- 
ating by OVERFLOW and coarsened by doubling the cell 
size for each successive level of grids that extended the far 
field computational domain. This procedure was used un- 
til the far field boundary was extended approximately 527? 
away from the hub center. The numerical simulations had 
a total of 160 million grid points, including 40 million grid 
points used for the body and rotor volume grids. 

Fuselage and Rotor Blade Geometries 

The ROBIN-mod7 geometry is shown in Figure 3 along 
with model dimensions for the joint U.S. Army and NASA 
forward flight test to be conducted in the 14-by 22-Foot 
Subsonic Tunnel at NASA Langley. The fuselage geometry 
for the forward flight test is the same as the isolated fuse- 
lage experiment of Schaeffler et al. (Ref. 10) with the added 
inclusion of an upper pylon that partially encloses the rotor 
shaft. The model for the isolated fuselage test had a fuse- 
lage length of 28.24” as compared to a length of 123.931” 
for the planned forward flight test fuselage. The Reynolds 
number, based on the fuselage length of the forward flight 
test, will range from 11-16 million, an order of magnitude 
larger than the 1.6 million of the previous isolated fuselage 
test. 


The rotor blades were acquired specifically for a PSP 
validation test and were used in a recent hover test of Wong 
et al. (Ref. 19) at the 14-by 22-Foot Subsonic Tunnel Rotor 
Test Cell. The blade uses Government RC-series airfoils 
with the planform shown in Figure 4. The rotor has a ra- 
dius of 66.5” and a linear twist of -14°. The blade has a 
chord length of 5.45” with a 30° tip sweep and a 3.27” tip 
chord length. Figure 4 show a notch in the planform from 
blade station 7.87 to 16.75. This notched out area contains 
a connector fairing to accommodate wiring from the blade 
pressure measurements on all four blades and is shown in 
Figure 4. The flap and lead-lag hinge are collocated 3.00” 
from the hub center. 

The ROBIN fuselage was originally developed using the 
radius, of the rotor paired with the fuselage as a parameter 
for the non-dimensional description of the fuselage geom- 
etry. The length of the fuselage is defined as twice the ro- 
tor radius and the cartesian coordinates used to define the 
fuselage surface are all non-dimensionalized by this paired- 
rotor radius. Over the years, the ROBIN fuselage has been 
paired in wind tunnel tests with rotors that are either under 
or oversized based upon this non-dimensionalization. The 
current research effort is not an exception to this. The PSP 
rotor is oversized for the current fuselage by 7.3%. The 
scaling radius for the fuselage, 7?^, is 61.9655 inches ver- 
sus the PSP rotor radius, 7?, of 66.50 inches. 

Flight Conditions 

The goal of this investigation is to reduce the fuselage drag 
at high speed cruise where the fuselage drag is a significant 
portion of the entire vehicle drag. An advance ratio of jli = 
0.35 was selected as a typical cruise condition for maximum 
range with Cj / (5 = 0.08. The hover tip Mach number for 
this investigation is 0.598 and the freestream Mach number 
is 0.2093. The shaft tilt angle of the rotor with respect to the 
fuselage is fixed at an angle of -3° for the wind tunnel model 
and hence was also fixed for all of the CFD simulations. An 
a =-3° was used as the primary simulation condition with 
some CFD simulations performed at a = -5°. The a = - 
3° and -5° corresponds to a rotor shaft angle, a s , of -6° 
and -8° respectively and are typical of rotorcraft cruise for 
maximum range. For comparison purposes, a = -3° and - 
5° for the isolated fuselage cases were also simulated with 
and without flow control in addition to an a = 0° baseline 
isolated fuselage case. 

Results 

Isolated Fuselage Baseline 

The fuselage without a rotor was simulated to better under- 
stand the effects of the rotor on the fuselage with and with- 
out flow control. The isolated fuselage simulations have a 
freestream Mach number of 0.2093 with a Reynolds num- 
ber based on the fuselage length of 14.9 million, matching 



the rotor/fuselage flow conditions. The isolated fuselage 
was evaluated at a = 0°, -3°, and -5° to compare with the 
rotor/fuselage simulations. The a = 0° case was simulated 
by solving the steady- state RANS equations and had good 
convergence in the RANS equation residuals. The steady - 
state solution for the a = -3° case did not converge as well. 
In this case, the drag had a sudden jump after many itera- 
tions to a different flow state and a much higher drag coef- 
ficient. This jump in the drag and lift convergence histories 
can be seen in Figs. 5 and 6 respectively. As seen in the 
drag history plot in Fig. 5, the drag coefficient converges to 
approximately 0.1 after 30,000 iterations. At this point the 
drag seems converged, however the lift history has not. Ad- 
ditional flow solver iterations result in a drag rise that even- 
tually converges after 75,000 iterations for both the drag 
and lift coefficients. 

This sudden drag rise between the flow states in the 
steady-state solution and the lack of convergence of the 
residual prompted a time-accurate flow analysis of the iso- 
lated fuselage at a = -3°. The unsteady flow simulation for 
the a = -3° case was initiated from a freestream condition 
with the time history of the drag and lift convergence pre- 
sented in Figs. 7 and 8. Initially the unsteady flow simula- 
tion was run using a non-dimensional time step of DT=0.58, 
where DT is non-dimensionalized by the hover tip velocity 
and the unit grid length of one inch. Using this time step, 
Cp converges to 0.1 as seen in the drag history plots. Look- 
ing at the unsteady simulation results, it was noticed that 
there was an unsteadiness at the juncture of the upper pylon 
and fuselage causing an unsteadiness in the juncture vortex 
that interacts with the flow in the ramp region. Decreasing 
DT by half to a value of 0.29, the unsteadiness of the junc- 
ture vortex was arrested while the flow in the ramp region 
became unsteady. This unsteady flow on the ramp increased 
the mean drag by 20% where the forces on the body become 
periodic after tUoo/2Rp = 30 as shown in Figs. 7 and 8. 

The isolated fuselage was then repositioned to a = -5° 
where the flow was initiated using the a = -3° flow solution. 
The results of the drag and lift histories are shown in Figs. 7 
and 8 starting at tUoo/2Rp = 40. The unsteady drag history 
indicates that the flow is fairly unstable and converges to a 
periodic flow state at tU 00 /2Rp = 67. The flow at this point 
shows the juncture vortex on the upper pylon to become un- 
steady when viewing a time lapse animation of the surface 
restricted streamlines on the fuselage. Cutting the time step 
in half to DT=0.145 did not affect the unsteadiness of the 
juncture vortex from the upper pylon or change the mean 
drag values. Based on these CFD results there seems to be 
a link between the unsteadiness of the juncture vortex and 
the unsteadiness on the ramp region that should be explored 
further. 

Surface streamlines and C p values for the baseline iso- 
lated fuselage are shown in Fig. 9 for the a = 0°, -3°, and 
-5° cases. One of the main differences between the three 
fuselage angles is the interaction of the upper pylon junc- 


ture vortex on the separated flow in the ramp region. In 
the previous research of Schaeffler et al. (Ref. 10), the fuse- 
lage was intentionally keep as clean as possible for the CFD 
validation effort and did not include the upper pylon. The 
juncture vortex for the a = 0° case has little interaction with 
the ramp flow separation as seen in Fig. 9(d, g, j). The a = 
-3° case in Fig. 9 does show the juncture vortex interact- 
ing with the ramp separation resulting in four different foci 
in the surface streamline pattern with a large suction pres- 
sure on the ramp edges. Overall, the ramp base pressure 
has decreased at a = -3°, as compared to the a = 0° case, 
contributing to the higher drag coefficient. Decreasing the 
fuselage angle-of-attack to a = -5°, resulted in the juncture 
vortex moving to a position higher on the ramp and increas- 
ing the base pressure 

Fuselage centerline C p plots for the isolate fuselage, for 
three different angles-of-attack, are given in Figs. 10 and 
1 1 for the lower and upper surface respectively. The cen- 
terline C p profile for the lower surface has a favorable pres- 
sure gradient on the nose with an adverse pressure gradi- 
ent starting at the geometric transition point from the nose 
to the straight body section. The pressure decreases as it 
approaches the ramp section with a suction peak just down- 
stream from x/Rp = 1.0. The flow separates at x/Rf = 1.07 
for a = 0° where the pressure flattens. The centerline pres- 
sures begin to increase atJv/T^/r = 1.15, where the ramp tran- 
sitions into the tail section. This increased pressure is most 
prominent for the -3° case and can also be seen in the ramp 
C p contours in Fig. 9(h). A close-up view of the pressure 
suction peak at the beginning of the ramp section is shown 
in Fig. 12 for the three a cases. This comparison indicates 
that there is only a small difference in the peak suction pres- 
sures for the different angles-of-attack with the a = -3° hav- 
ing the largest suction pressure. The upper centerlines C p 
values in Fig. 1 1 have a similar favorable pressure gradient 
on the nose section with a large adverse pressure gradient 
due to the upper pylon. Overall, there is little difference 
between the three a cases for the upper and lower surfaces 
with most of the differences on the nose section and where 
the ramp transitions into the tail. 

Isolated Fuselage with Flow Control 

A series of eight slots that form a nearly continuous U- 
shape on the ramp region of the fuselage was found to 
perform best out of the other slot configurations evalu- 
ated in the isolated fuselage research by NASA and ON- 
ERA (Ref. 10). This slot configuration is shown in Fig. 3 
and was generated by locating the slots on the intersection 
of a plane at a constant x/Rp = 1.155. The intersecting 
plane was rotated 23° at x/Rp = 1.155 and z/Rf = 0 along 
the span wise axis as shown in Fig. 3. This slot layout re- 
sulted in slots that are slightly downstream of the flow sep- 
aration line on the fuselage ramp. Building on the prior 
isolated fuselage research, we will now evaluate the perfor- 



mance of the flow control for the fuselage and upper pylon 
with and without the rotor. The primary slot height for this 
investigation is 0.0195”. The jet is set to an angle of 25° to 
the local surface tangent and blows downstream and normal 
to the local slot span wise vector. This jet angle was found 
to perform well by Schaeffler et al. (Ref. 10). The previous 
research did show a slight increase in performance for a 15° 
jet slot angle; however, this slot angle was difficult to im- 
plement experimentally. A larger slot height of 0.030” was 
also simulated for a select number of cases as guidance for 
the design of the wind tunnel model flow control slots. 

The flow control simulations were initiated from their 
respective baseline cases using ZNMF jets with a non- 
dimensional frequency of F + =0.56. This frequency cor- 
responds to four jet blowing and suction cycles per one ro- 
tor revolution in the rotor/fuselage simulations discussed 
below. Based on the results from the NASA and ON- 
ERA (Ref. 10) research, each blowing slot is 180° out-of- 
phase from its adjacent slots. This was found to greatly 
minimize oscillations in the fuselage drag and provide the 
best drag reduction performance. 

a = -3° Case The isolated fuselage drag history results, 
using the 0.030” slots for a VR of 1.5, resulted in a mean Cd 
value of 0.094 or a 22% reduction from the baseline drag. 
This case has a peak C p value of 0.0068, which is based 
on the maximum cross-sectional area, Acs , of the fuselage 
body, excluding the upper pylon. The peak C p value is cal- 
culated using the total mass flux exiting all the slots and the 
peak exit velocity. Since half the jets are 180° out-of-phase 
from each other, only half the jets are blowing during the 
peak condition while the other half are at peak suction. The 
smaller 0.0195” slot height, for the same VR of 1.5, has 
a C p value 0.0045 and results in an overall mean drag re- 
duction; however, this drag history switches intermediately 
between a Cd value of 0.10 and 0.1 1 as seen in Fig. 13. In- 
creasing the VR to 2.0 for the 0.0195” slot, increases C p to 
0.008 and removes the large unsteadiness in the drag his- 
tory. The higher VR of 2.0 results in a mean drag value of 
0.099, which is a 17.5% reduction from the baseline drag. It 
is noted that the larger slot, which has a C p value of 0.0068, 
performed better than the smaller slot at the higher C p value 
of 0.0080. Similar results were seen in the CFD simulations 
for the low Reynolds number investigations (Ref. 10) for 
steady blowing. Additional flow control simulations with 
the larger slot height needs to be performed in order to bet- 
ter understand the performance difference between the two 
slot sizes for the isolated fuselage. Discussions of the per- 
formance between the two slot sizes are made below in the 
Rotor/Fuselage with Flow Control section. 

The lift results for the a = -3° isolated fuselage case in 
Fig. 14 show a reduction of the download on the fuselage. 
The larger 0.030” slot height with a VR of 1.5 has a mean 
Cl value of -0.088 as compared to the baseline mean Cl 
value of -0.12 and is a 28% reduction of the fuselage down 


load. The smaller 0.0195” slot for a VR of 1.5 had a similar 
mean value; however, the lift time history shows the same 
intermittent flow switching seen in the drag. Increasing the 
VR to 2.0, for the 0.0195” slot, results in a mean Cl value of 
-0.114 ora7.3% reduction of fuselage baseline download 
force. 

To better understand the effects of the flow control on the 
fuselage drag, a view of the surface restricted streamlines 
and C p contours for the baseline and flow control cases are 
shown in Fig. 15. The flow control case in this figure has a 
slot height of 0.0195” and a VR of 2.0. This figure shows 
a snapshot of the flow on the ramp for a jet phase of (]) = 
0°, 90°, and 270°. At (f) = 0°, seen in Fig. 15(b), all the 
jets have zero velocity. Comparing the C p contours to the 
baseline case in Fig. 15(a), it can be seen that the ramp base 
pressure has significantly increased, contributing to a fuse- 
lage drag reduction. This figure also shows a reduction of 
the flow separation on the ramp with an increase of the suc- 
tion pressure on edges of the ramp region. This increase 
in the suction pressure is due to the flow being reattached 
on the ramp, increasing the external flow velocity on the 
ramp edges. This suction pressure does contribute to in- 
creasing the fuselage drag, reducing the effectiveness of the 
flow control jets. At (|) = 90° and 270°, seen in Fig. 15(c) 
and (d), half the jets are at peak blowing and the other half 
at peak suction. These two jet phases are fairly similar to 
Fig. 15(b) in that the ramp base pressure is higher than the 
baseline case and the flow separation locations have move 
inward with a large suction pressure on the comers of the 
ramp. These two figures also show a high pressure region 
just downstream of the slots that are at peak blowing. 


a = -5° Case Decreasing the fuselage a to -5° resulted 
in a loss of performance for the ZNMF flow control. Only 
the 0.0195” slots were simulated for this case with the drag 
results shown in Fig. 16. The VR of 1.86 and 2.5 both in- 
creased the mean drag coefficient with a marginal drag re- 
duction for a VR of 3.0. While the ramp base pressure has 
improved for these cases, the suction pressure on the cor- 
ners of the ramp have also increased. From these results it 
is apparent that the performance of the flow control is af- 
fected by the rotation of the fuselage to a more nose down- 
ward position. This rotation results in the suction pressure 
on the bottom of the ramp having a larger component in the 
drag direction, reducing the effectiveness of the flow con- 
trol. The loss of performance for the fuselage lift is also 
seen in Fig. 17 where the download force on the fuselage 
has doubled. The nose down rotation of the fuselage also 
means that the increase base pressure on the ramp also has 
an increased downward component, contributing to the total 
download on the fuselage. 

The surface C p contours and streamlines for the ramp re- 
gion, are shown in Fig. 18 for the a = -5° isolated fuselage 
case. Like the a = -3° case, the flow control is showing a 



similar increased ramp base pressure with an increased suc- 
tion pressure on the corners on the rear ramp. A comparison 
of the C p contours between the a = -3° and -5° cases indi- 
cates similar performance for the flow control in terms of 
increasing the base pressure on the ramp. 

To better understand the increased fuselage download 
using AFC for the a = -5° case, a contour plot of the differ- 
ence between the surface C p for the AFC and the baseline 
case, AC P , are made in Figs. 19 and 20. These two figures 
show a bottom view of the fuselage for a = -3° and -5°. 
Comparing the two contour plots reveals that the a =-5° 
AFC case has a larger increase in the suction pressure near 
the ramp comers and a significant region of increased lower 
pressure on the tail boom, as compared to a = -3°. There- 
fore, the a = -5° contour plot identifies the source of the 
download increase as the low pressure region on the lower 
ramp corners and tail boom. 

Isolated Fuselage Centerline Pressures The isolated 
fuselage centerline C p profiles on the lower surface, with 
and without flow control, are shown in Figs. 21 and 22. 
These two figures show how the flow control increases the 
pressure on the centerline of the ramp as compared to the 
baseline case. These figures also show an increase of the 
suction pressure on the edge of the ramp region. Note that 
the flow control has also resulted in a lower pressure on the 
tail boom end cap with the a = -5° case having a significant 
decrease in pressure. 

Rotor/Fuselage Baseline 

The simulation of the rotor and fuselage was performed at 
the forward fight conditions, described in the Flight Con- 
ditions section above, for an advance ratio of ju = 0.35 and 
rotor thrust of CT /a = 0.08, at a = -3° and -5°. The rotor 
dynamics were simulated using a loose coupling between 
the flow solver and the comprehensive code, CAMRAD-II 
(Refs. 20,21). The rotor blades are modeled as rigid blades 
with flapping and lead/lag motions. In the rotor/fuselage 
simulations, a target rotor thrust is set and the comprehen- 
sive code adjusts the collective and cyclic inputs to match 
the commanded thrust while minimizing flapping. The 
blade motions from the comprehensive code are then used 
in the unsteady CFD simulations. The aerodynamic forces 
on the rotor blades are then passed back to the comprehen- 
sive code and used to update the dynamic motion of the 
rotor blades and to trim the rotor, matching the target thrust 
while minimizing flapping. This coupling was performed 
every half revolution for the rotor/fuselage simulations. 

The baseline fuselage drag and lift for the rotor fuselage 
simulations are shown in Figs. 23 and 24 for both a = -3° 
and -5°. The a =-3° baseline case was simulated for 9 ro- 
tor revolutions and the a =-5° case for 16 rotor revolutions 
allowing the fuselage drag and lift to converge. The time 


history showed the fuselage lift took longer than the drag to 
reach a converged condition. For both cases, the rotor lift 
converged within 4 revolutions. The unsteady fuselage drag 
has large 4/rev oscillations that are attributed to a pressure 
increase on the front of the fuselage when a rotor blades 
passes directly over the fuselage. A comparison of the fuse- 
lage drag between the a = -3° and -5° cases shows an in- 
crease in the amplitude of the drag for a = -5°. However, 
the mean drag coefficient for both a cases was 0.18, a 50% 
drag increase over the isolated fuselage mean Cd value of 
0.12. The mean lift coefficient for the a = -3° case was 
-0.24, a 50% increase over the isolated fuselage value of 
-0.12. The lift for a = -5° in Fig. 24 is slowly increasing 
due to a low frequency oscillation seen in the fuselage lift. 

Surface restricted streamlines and C p contours for the 
baseline rotor/fuselage simulations are shown in Figs. 25 
and 26 at three different blade azimuth angles for a =-3° 
and -5° respectively. These figures show that the flow sep- 
aration, with respect to the slot locations, does not vary 
greatly from the baseline isolated fuselage cases. This indi- 
cates that the flow control slots do not need to be relocated 
for the rotor/fuselage case. 

The baseline fuselage centerline C p profiles are shown 
in Figs. 27-30 for the a = -3° and -5° rotor/fuselage cases. 
These figures show the C p profiles for three blade azimuth 
angles of \|/ = 0°, 30°, and 60°. The lower surface C p pro- 
files for both cases show little difference for the three rotor 
azimuth angles. The upper surface C p profiles are similar 
between the \\f = 30° and 60°, with the \|/ = 0° case hav- 
ing a higher pressure on the nose and tail as a result of the 
blade being directly over these regions. Therefore, most of 
the unsteady pressure forces are related to the passage of 
the rotor blades over the nose and tail section as expected. 
Comparing the lower surface C p profiles between the a = 
-3° and -5° cases, shows the a = -3° case to have a larger 
suction peak on the ramp section with a higher pressure on 
the upper part of the ramp as compared to a = -5°. This 
higher pressure was also seen in the C p surface contours in 
Figs. 25 and 26. 

Rotor/Fuselage with Flow Control 

Simulation of the rotor/fuselage with flow control was per- 
formed using the same slot location as in the isolated fuse- 
lage case. The ZNMF jets had a non-dimensional frequency 
of F+ =0.56, resulting in 4 blowing and suction cycles per 
rotor revolution. Slot heights of 0.0195” and 0.030” were 
evaluated for a = -3° with varying VR values, where the 
a = -5° was only evaluated using the 0.0195” slot height. 

The fuselage unsteady drag and lift coefficients, as a 
function of the rotor azimuth angle, are given in Figs. 31- 
34 for a = -3° and -5°. These comparisons show that the 
flow control resulted in an overall drag decrease, shifting 
the baseline unsteady drag profile downward. Likewise, the 



flow control resulted in an overall increase in the fuselage 
lift coefficient for a =-3° and a decrease for a = -5°. 

The flow control performance of all of the CFD simula- 
tions are presented in Figs. 35-38. The flow control perfor- 
mance for the two jet slot heights at a = -3° are shown in 
Fig. 35. These results show that the two slot heights for the 
same C ' p value, only have a small difference with the larger 
0.030” slot having slightly improved performance. A com- 
parison of the two slots for the fuselage lift in Fig. 36, show 
the larger slot performing approximately 5% better than the 
smaller slot. Therefore, the flow control performance for 
the fuselage drag shows an expected dependence on C p , 
while the lift indicates a dependence on the slot height and 

c r 

The percentage change of Cp and Cl from the baseline 
cases for increasing C p are presented in Figs. 37 and 38. 
The percentage mean drag change in Fig. 37 indicates that 
the isolated fuselage and rotor/fuselage for a = -3°, have 
a similar percentage drag reduction from their respective 
baseline values. The a = -5° performance results for the 
isolated fuselage and the rotor/fuselage cases are not similar 
as was seen for the a = -3° case. 

The fuselage lift performance in Fig. 38 indicates that 
both the isolated fuselage and rotor/fuselage simulations 
initially show an increase in lift on the fuselage for increas- 
ing C p at the a = -3° condition. The isolated fuselage lift 
shows a decrease as C p increases from 0.0045 to 0.0080, 
whereas the lift for the rotor/fuselage simulation continues 
to increase. The a = -5° fuselage flow control case for the 
rotor/fuselage simulation has an initial decrease of fuselage 
lift as C/j increases. The lift performance then recovers as 
C p is increased further from 0.0045 to 0.0120, approach- 
ing the original baseline lift value. The isolated fuselage 
lift results are significantly different from the rotor/fuselage 
simulations at a =-5° showing a much larger percentage 
download increase from the baseline for increasing C p . 

The effect of the flow control on the ramp surface pres- 
sure can be seen in Figs. 39 and 40 for the 0.0195” slot 
height and VR = 1.86. These two figures show the ramp 
region for three different blade azimuth angles. Overall, 
the surface C p contours show an increase of the ramp base 
pressure with AFC, reducing fuselage drag. The flow sepa- 
ration locations have also moved inward with a larger suc- 
tion pressure on the edges of the ramp region, contributing 
to an increase in fuselage drag. 

A comparison of the spanwise surface C p values halfway 
up the ramp, on a constant z/Rf = -0.075 plane, are made in 
Figs. 41 and 42 for both the a = -3° and -5° rotor/fuselage 
cases. The instantaneous C p plots are given for \|/ = 0°, 30°, 
and 60° . A comparison of the baseline C p values for both 
cases reveals similar pressure profiles for the three blade az- 
imuth angles where the a = -5° baseline case has a slightly 
higher base pressure than the a = -3° baseline case. The C p 
profiles with AFC, for both a cases, have an increase in the 


ramp base pressure with an increase of the suction pressure 
on the corners. 

Figures 43 and 44 compare the lower surface centerline 
C p profiles for the baseline and flow control cases for the 
two angle-of-attack cases at an azimuth angle of \\f = 0°, 
which is characteristic of the other blade azimuth angles. 
The upper C p centerline profile was not effected by the flow 
control, however flow control did modify the lower center- 
line C p profile near the ramp region. Both a cases show an 
increase in the suction pressure and higher ramp base pres- 
sure with AFC. 

A view of the fuselage and rotor wakes are presented 
in Figs. 45-48 for the baseline and flow controlled CFD 
simulations showing iso-surface contours of the Q-criterion 
shaded by C p . Figures 45 and 46 show the rotor tip vor- 
tices and the juncture vortex from the upper pylon. Fig- 
ures 47 and 48 are side views of the simulation, highlighting 
the interaction of the rotor wake vorticity with the fuselage. 
This view indicates that there is little interaction of the ro- 
tor wake vorticity with fuselage ramp region at this forward 
flight condition. It also shows the trajectory of the juncture 
vortex on the side of the fuselage, entering the ramp re- 
gion. The side view highlights the modification of the wake 
from the separated flow on the ramp and unsteadiness of the 
streamwise vortex from the ramp region. 

Summary 

This investigation explores the effectiveness of flow con- 
trol for rotorcraft fuselage drag reduction using ZNMF jets. 
This research extends the isolated fuselage work of Scha- 
effler et al. (Ref. 10) to a more relevant environment by in- 
cluding a rotor. This effort is also complementary to a future 
1/3 scaled rotor/fuselage wind tunnel experiment planned at 
the NASA 14-by 22-Foot Subsonic Tunnel. The outcome 
of this numerical investigation will be used as guidance for 
the design of a fuselage flow controlled model for the future 
test. 

The isolated fuselage simulation used the same fuselage 
geometry as in the previous investigation of Schaeffler et 
al. (Ref. 10) with the addition of an upper pylon that par- 
tially encloses the rotor shaft. It was discovered that the 
inclusion of the upper pylon at a = -3° resulted in a junc- 
ture vortex that interacted with the separated flow in the 
ramp region. This interaction of the juncture vortex with the 
ramp resulted in unsteady vortex shedding from the ramp, 
increasing the fuselage drag by 20%. The isolated fuselage 
was simulated at a typical long range cruise condition were 
Moo = 0.209 with the fuselage at a = -3°. A RANS sim- 
ulation of the isolated fuselage with flow control resulted 
in a 17% drag reduction from the baseline with a 7% de- 
crease of the fuselage download force. The flow control at 
this condition had a slot height of 0.0195” and a peak jet 
VR of 2.0, resulting in a C p value of 0.008. At a = -5°, 
the flow control resulted in a drag increase. However, a 



marginal drag reduction of 2% was seen for a very large C ^ 
value of 0.015. In addition to the drag increase, the fuse- 
lage download also increased dramatically for the -5° flow 
control case by 140% for C ^ = 0.0069. The flow control for 
the a = -5° case resulted in lower pressures on the bottom 
ramp comer and the bottom of the tail boom, as compared 
to the baseline pressure. These low pressure regions were 
identified as the source of the increase fuselage download. 

Including the rotor in the numerical simulations intro- 
duced a large 4/rev oscillation of the fuselage forces, in- 
creasing the mean fuselage drag by 50% with a 50% in- 
crease of the fuselage download as compared to the iso- 
lated fuselage at a = -3°. Turning on the fuselage flow con- 
trol, for the 0.0195” slot with a jet VR of 1.86 (C AI = 0.0069), 
resulted in a 17% drag reduction with a 22% reduction of 
the fuselage download force. The reduction of the fuselage 
download has a compounding positive effect on the over- 
all vehicle performance, since the rotor thrust can be re- 
duced, which further reduces the download on the fuselage. 
Evaluation of two different flow control jet slot heights of 
0.0195” and 0.030” resulted in similar drag performance 
curves when plotted as a function of C^. These result indi- 
cated that the performance of the flow control for this con- 
figuration is a function of the jet momentum coefficient, C^. 

Comparing the flow control performance between the 
isolated fuselage and rotor/fuselage simulations at a = -3° 
showed that the drag performance, in terms of percentage 
drag reduction, was similar over the range of C ^ evaluated. 
Unfortunately, this was not the case for a = -5°, where the 
isolated fuselage, unlike the rotor/fuselage case, had an ini- 
tial drag increase with increasing flow control and a sig- 
nificantly larger increase of the fuselage download. 
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Fig. 1. Surface grids on the fuselage and rotor. 


Fig. 2. Detail of grids used for flow control slots. 



Fig. 3. Fuselage geometry and slot locations used for the CFD simulations where Rf is 61.9655 for the 1/3 scale 14- 
by 22-foot Subsonic Wind Tunnel model 



Fig. 4. Blade layout. 








Fig. 5. Isolated fuselage baseline drag coefficient con- Fig. 6. Isolated fuselage baseline lift coefficient conver- 
vergence history for steady-state solution. gence history for steady-state solution. 




Fig. 7. Drag coefficient history for the time-accurate Fig. 8. Lift coefficient history for the time-accurate sim- 

simulation of the baseline isolated fuselage for a = -3° ulation of the baseline isolated fuselage for a = -3° and 

and -5°. -5°. 



Fig. 9. Baseline fuselage for Moo=0.2093 at a = 0°, -3°, and -5° shaded by the surface C p with surface restricted 
streamlines and flow control slots outlined in red: (a-c) side view, (d-f) view showing juncture vortex from upper 
pylon interacting with ramp (g-i) rear view of ramp region (j-1) close-up of the side view. 




Fig. 10. Fuselage centerline C p values for the isolated Fig. 11. Fuselage centerline C p values for the isolated 
fuselage lower surface at a = 0°, -3°, and -5°. fuselage upper surface at a = 0°, -3°, and -5°. 



Fig. 12. Fuselage centerline C p values for the isolated fuselage lower surface at a = 0°, -3°, and -5 



Fig. 13. Drag coefficient history for the time-accurate 
simulation of the isolated fuselage with ZNMF control 
jets for a = -3°. All simulations are using a jet blowing 
frequency of F + = 0.56. 
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Fig. 14. Lift coefficient history for the time-accurate 
simulation of the isolated fuselage with ZNMF control 
jets for a = -3°. All simulations are using a jet blowing 
frequency of F + = 0.56. 
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Fig. 15. Comparison of the C p contours and surface restricted streamlines on the ramp region for the isolated 
fuselage at a = -3° for (a) the baseline case (i.e. no flow control) and AFC case at three actuator phases, (b) § = 0°, 
(c) <|> = 90°, and (d) 4> = 270° for a VR = 2.0. 
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Fig. 16. Drag coefficient history for the time-accurate 
simulation of the isolated fuselage with ZNMF control 
jets for a = -5°. All simulations are using a jet blowing 
frequency of F + = 0.56. 


Fig. 17. Lift coefficient history for the time-accurate 
simulation of the isolated fuselage with ZNMF control 
jets for a = -5°. All simulations are using a jet blowing 
frequency of F + = 0.56. 




Fig. 18. Comparison of the C p contours and surface restricted streamlines on the ramp region for the isolated 
fuselage at a = -5° for (a) the baseline case (i.e. no flow control) and AFC case at three actuator phases, (b) § = 0°, 
(c) 4> = 90°, and (d) 4> = 270° for a VR = 2.5. 
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Fig. 19. Bottom view of the computed C p difference between the baseline and AFC isolated fuselage at a = -3° case. 
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Fig. 20. Bottom view of the computed C p difference between the baseline and AFC isolated fuselage at a = -5° case. 
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Fig. 21. Fuselage centerline C p values for the isolated 
fuselage lower surface at a = -3° for the baseline and 
flow control cases. 
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Fig. 22. Fuselage centerline C p values for the isolated 
fuselage upper surface at a = -5° for the baseline and 
flow control cases. 
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Fig. 23. Baseline fuselage drag coefficient for the ro- 
tor/fuselage simulation at a = -3° and -5° with M «, = 

0.209, and CT/a = 0.08. 
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Fig. 24. Baseline fuselage lift coefficient for the ro- 
tor/fuselage simulation at a = -3° and -5° with M «, = 

0.209, and CT/a = 0.08. 
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( C ) C p : - 0.8 - 0.6 - 0.4 - 0.2 0 0.2 0.4 

Fig. 25. Rear ramp view of the baseline rotor/fuselage 
simulation for the a = -3° case. Surface C p contours 
with surface restricted streamlines at (a) \\f = 0°, (b) \|/ = 
30°, and (c) \| / = 60°. (Simulation parameters: ju = 0.35, 
CT /a = 0.08, Afoo= 0.2093) 
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Fig. 26. Rear ramp view of the baseline rotor/fuselage 
simulation for the a = -5° case. Surface C p contours 
with surface restricted streamlines at (a) \|/ = 0°, (b) \j/ = 
30°, and (c) \]/ = 60°. (Simulation parameters: \i = 0.35, 
CT /a = 0.08, Afoo= 0.2093) 





Fig. 27. Lower fuselage centerline C p values for the 
baseline rotor/fuselage simulation at a = -3° for vary- 
ing rotor azimuth angles. 
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Fig. 29. Lower fuselage centerline C p values for the 
baseline rotor/fuselage simulation at a = -5° for vary- 
ing rotor azimuth angles. 



Fig. 28. Upper fuselage centerline C p values for the 
baseline rotor/fuselage simulation at a = -3° for vary- 
ing rotor azimuth angles. 
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Fig. 30. Upper fuselage centerline C p values for the 
baseline rotor/fuselage simulation at a = -5° for vary- 
ing rotor azimuth angles. 



Fig. 31. Unsteady fuselage drag coefficient of the ro- 
tor/fuselage simulations at a = -3° for varying jet VR 
values. (Simulation parameters: F + = 0.56, ju = 0.35, 
CT/a = 0.08, Moo=0.2093) 



Fig. 32. Unsteady fuselage lift coefficient of the ro- 
tor/fuselage simulations at a = -3° for varying jet VR 
values. (Simulation parameters: F + = 0.56, ju = 0.35, 
CT/a = 0.08, Moo=0.2093) 
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Fig. 33. Unsteady fuselage drag coefficient of the ro- 
tor/fuselage simulations at a = -5° for varying jet VR 
values. (Simulation parameters: F + = 0.56, ju = 0.35, 
CT/a = 0.08, Moo=0.2093) 
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Fig. 34. Unsteady fuselage lift coefficient of the ro- 
tor/fuselage simulations at a = -5° for varying jet VR 
values. (Simulation parameters: F + = 0.56, ju = 0.35, 
CT/a = 0.08, Moo=0.2093) 


0.00 



peak 


U 

c 

(0 

QJ 


-0.05 


- 0.10 


-0.15 


- 0.20 


-0.25 


-0.30 




-A- Rotor/Body; -3 deg. AoA; 0.030" slot 
-□-Rotor/Body; -3 deg. AoA; 0.0195" slot 
-0- Rotor/Body; -5 deg. AoA; 0.0195" slot 
-O-Body; -3 deg. AoA; 0.0195" slot 
-A-Body; -5 deg. AoA; 0.0195" slot 


















A 










0.000 0.005 0.010 0.015 0.020 

C, peak 


Fig. 35. Mean fuselage drag coefficient as a function 
of C,j for the isolated fuselage and rotor/fuselage CFD 
simulations. 


Fig. 36. Mean fuselage lift coefficient as a function of C ^ 
for the isolated fuselage and rotor/fuselage CFD simula- 
tions 



0.000 0.005 0.010 0.015 0.020 

C, peak 



peak 


Fig. 37. Percentage change in the mean fuselage drag Fig. 38. Percentage change in the mean fuselage lift co- 
coefficient as a function of for the isolated fuselage efficient as a function of for the isolated fuselage and 

and rotor/fuselage CFD simulations. rotor/fuselage CFD simulations. 
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Fig. 39. Rear ramp view of the baseline rotor/fuselage 
simulation for a = -3°. Surface C p contours with sur- 
face restricted streamlines at (a) \\f = 0°, (b) \\f = 30° , 
and (c) \| / = 60°. (Simulation parameters: VR=1.86, h = 
0.0195”, C p = 0.0069, ju = 0.35, CT /o = 0.08, 0.2093) 
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( C ) C p : - 0.8 - 0.6 - 0.4 - 0.2 0 0.2 0.4 

Fig. 40. Rear ramp view of the baseline rotor/fuselage 
simulation for a = -5°. Surface C p contours with sur- 
face restricted streamlines at (a) \\f = 0°, (b) \\f = 30°, 
and (c) \| f = 60°. (Simulation parameters: VR=1.86, h = 
0.0195”, C p = 0.0069, ju = 0.35, CT /a = 0.08, M«f= 0.2093) 
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Fig. 41. Spanwise surface C p on the ramp for the base- 
line and AFC cases at a = -3°. 
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Fig. 43. Fuselage centerline C p values for the ro- 
tor/fuselage simulation at a = -3° comparing the base- 
line and AFC cases. 
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Fig. 42. Spanwise surface C p on the ramp for the base- 
line and AFC cases at a =-5°. 
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Fig. 44. Fuselage lower centerline C p values for the ro- 
tor/fuselage simulation at a = -5° comparing the base- 
line and AFC cases. 



Fig. 45. Baseline (no flow control) case showing an iso- 
surface contour of the Q-criterion shaded by C p for 
forward flight for a =-3°, £/«, = 0.2093,^ = 0.35, and 
C T /o = 0.08. 



Fig. 47. Baseline (no flow control) case showing a side 
view of an iso-surface contour of the Q-criterion shaded 
by Cp for forward flight at a =-3°, ju = 0.35, Uoo = 
0.2093, and C T /o = 0.08. 



Fig. 46. AFC case showing an iso-surface contour of the 
Q-criterion shaded by C p for forward flight for a =-3°, 
ju = 0.35, Uoo = 0.2093, and C r /o = 0.08. 



C p : - 0.5 - 0.4 - 0.3 - 0.2 - 0.1 0 0.1 0.2 0.3 0.4 0.5 


Fig. 48. AFC case showing a side view of an iso-surface 
contour of the Q-criterion shaded by Cp for forward 
flight at a =-3°, ju = 0.35, Uoo = 0.2093, and C T /a = 0.08. 


